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Abstract 

In the lattice work by Miller [H [2] and in the work by Zwanziger [3] a 
linear growth of the trace anomaly for high temperatures was found in pure 
SU(2) and SU(3) Yang-Mills theories. These results show the remarkable 
property that the corresponding systems are strong interacting even at high 
temperatures. We show that within an analytical approach to Yang- Mills 
thermodynamics this linear rise is obtained and is directly connected to the 
presence of a temperature-dependent ground state, which describes (part of) 
the nonperturbative nature of the Yang-Mills system. Our predictions are in 
approximate agreement with [H EJ [3] . 



1 Introduction 



Studies of pure SU(2) and SU(3) Yang-Mills theories at finite temperature T repre- 
sent an important subject in modern high energy physics and a fundamental step 
toward the analysis of the properties of the Quark-Gluon-Plasma. In particular, the 
question if and for which temperatures a gas of (almost) free gluons offers a valid 
approximate description of such strong interacting theories goes to the heart of the 
problem. A key-quantity to address this issue is the trace of the stress energy tensor 
= p — 3p, where p and p are the energy density and the pressure of the system. In 
fact, the trace 9^ vanishes for a free gas of massless gluons, being a consequence of 
conformal invariance of Yang-Mills theories. Thus, one would expect that 9^ van- 
ishes in the high-temeperature limit, where a Yang- Mills theory is asymptotically 
free. 

In the lattice simulations of Refs. [TJ [2j H], based on the lattice works of Refs. 
[5], a striking result was found: the trace of the stress-energy tensor 9 m = p — 3p 
of pure SU(2) and SU(3) Yang-Mills thermodynamics in the deconfining phase does 
not vanish (and neither decreases) but grows linearly with T for T > 2 T c where T c 
denotes the phase boundary. As emphasized in [TJ , this result means that the system 
is strong interacting even at high T and puts strong constraints on models aiming at 
describing these Yang-Mills theories at finite T. Notice that in [6] a quadratic rise 
at high T (instead of linear) of 9^ is inferred out of lattice data in the SU(3) case. 
On one hand this difference points to the difficulty in extracting high T behavior 
for thermodynamical quantities, on the other hand it confirms that nonperturative 
physics is at work in the high temperature domain. 

While phenomeno logical models [TJ [H [9] predict that, in accord with asymptotic 
freedom [ID] , the interaction measure A = vanishes for T — > oo, there is no 
general agreement on how fast A approaches zero. Namely, as discussed above, not 
only A but the quantity p—3p itself should vanish for T — > oo if the thermodynamical 
system is described by a noninteracting realtivistic gas of gluons. Notice that in the 
framework of resummed perturbation theory nonzero contributions to 9^ start at 
order g A [11], where it is shown that 9^ keeps growing even for high T. 

Assuming the correctness of the results of pQ, it is therefore important to study 
models which do describe the linear growth of 9^ in order to understand the ori- 
gin of this property. In this work we show that within the analytical approach to 
Yang-Mills thermodynamics of Ref. p2] the linear growth of p — 3p is predicted with 
slopes compatible with lattice results. This property is a consequence of the presence 
of a thermal ground state, which takes into account the contributions of calorons 
(i.e. instantons at finite T). Intuitively, this thermal ground state, described by 
an adjoint- scalar field 0, represents an average over such topologically non-trivial 
configurations. From a phenomenological perspective the ground-state contribution 
to p behaves like a temperature dependent bag constant, see for example [TJ. Fur- 
thermore, the introduced field <fi generates a Higgs mechanism: (some of) the gauge 
modes acquire a temperature-dependent mass, see details in following Section. 
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The ground-state contribution to the pressure p gs is negative (like the bag pres- 
sure is) and dominates when T \ T c . Lattice simulations of p do not agree on this 
point. Namely, within the differential method, close to T c negative pressure is pre- 
dicted [13] while the integral method [14] is designed to generate positive pressure. 
However, in the present report we investigate W only for T > 1.5 T c , thus well 
above the phase boundary. 

In Ref. [3] a linear growth of 9^ was obtained based on the suppression of the 
infrared modes by invoking a modified dispersion relation of Gribov type. The latter 
circumvents the well-known problems inherent in a perturbative treatment [T5j . In 
the alternative approach of [121 EE EL71 EH] the stability of the infrared sector is 
guaranteed by the above mentioned thermal ground state, embodied by an adjoint 
scalar field 0, taking into account nontrivial-topology fluctuations. Remarkably, 
this approach predicts the slopes a in m = aT once the Yang-Mills scale A (or 
alternatively T c ) is fixed. 

The present analysis is not intended to discuss phenomenological fits to lattice 
curves (for work in that respect see [3 O [8] and the review [19J). Rather, our goal 
is to show that the (asymptotically) linear growth of MAt is related to the presence 
of a nontrivial and temperature-dependent ground-state [121 HE] . 

The paper is organized as follows: First we briefly review the approach of Ref. [12] 
and list expressions for p and p in the SU(2) case. Subsequently, we show numerically 
and analytically that within this approach 0^ = aT for T ^> T c and compare a 
with the results of [H El [3]. We then discuss the SU(3) case and give our conclusions. 

2 The SU(2) case 

In Ref. [12] an adjoint scalar field <ft is derived which, upon spatial coarse graining, 
describes the dynamics of nontrivial, BPS-saturated topological configurations [201 
|2~T] giving rise to an energy density and pressure of the thermal ground state. The 
field acts as a background to the dynamics of topologically trivial fluctuations. 

The modulus \<f>\ = -y/^pf exhibits a dependence on T and on A. The cor- 
responding potential is V^(|0| 2 ) = A 6 /|0| 2 [12], [16] - In unitary gauge the effective 
Lagrangian describing a thermalized, pure SU(2) Yang-Mills system reads: 

CT = C [a,] = \ (GrW) 2 + 2e 2 |0| 2 ((ay 2 + (a 2 ) 2 ) + 2^ , (1) 

where e = e(T) denotes the temperature dependent effective gauge coupling (see 
below) which enters both into the effective field strength G a ^ iv and into the mass m 
for the fields a}; 2 . One has 

m 2 = m (T) 2 = ml = m\ = 4e 2 \(f)\ 2 , m\ = . (2) 

that is a temperature-dependent mass is generated in a nonperturbative fashion by 
the adjoint- scalar field 0. From the effective Lagrangian ([1]) one derives, on the one- 
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loop level (accurate to about 0.1% [T8J [22]), the energy density p and the pressure 
p as 

P = P3 + Pl,2 + Pgs, P = P3+ Pl,2 + Pgs , (3) 

where p\^ denotes the sum over the two massive modes. Explicitly we have: 



7T 



2 



dk k 2 Vm 2 + k 2 „ A 6 



vr 2 m4 f 00 dkk 2 , / 



p 3 = 2— T , p 1>2 = -6 T ^ — In ^1 - e , Pgs = -p gs . (5) 

Let us now rewrite the system in terms of dimensionless quantities: 

p=jt, p=f I , A = — a(A) = ^ = 2A_i |0| , (6) 

where the function a = a(A) is introduced for later use. In dependence of the 
dimensionless temperature A the dimensionless energy density p and pressure p read: 



_ vr 2 _ 3 /"» gyg+g _ 2(2vr) 4 

2 q poo 

p 3 = 2-, p^ = --j dxxHnil-e-^*) , p gs = - p gs . (8) 
By virtue of the Legendre transformation 

the evolution equation for a = a(A) follows: 

6A 3 A da \ „ . . 

1 = -(2^( A dA +a ) aBW - (10) 
Z" 00 1 

There exits a low-temperature attractor to the evolution described by Eq. (11 01) with 
a logarithmic pole at the critical temperature A c = 13.89 [12]. At A c the fluctuations 
a*' 2 decouple thermodynamically. The effective coupling is given as e = e(A) = 
a(A)A 3 / 2 /47r. For a C 1 the coupling e is constant: e = \/8tt. In fact 

"(A) = ^ (12) 

then is a solution to Eq. (JT0|) . This plateau is reached rapidly for increasing A, see 
figures in [12] . 
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Figure 1: The quantities ^ (left panel, gray curve: large-T behavior; black curve: 
finite-T behavior) and h (right panel) as functions of j- = |r in the case of SU(2). 



3 Linear growth of 9^ 

We now turn to the linear rise of 9^ within the approach described in the previous 
Section. In therms of dimensionless quantities one has 



&un P ~ 3p ( A x ' 



A 4 A 4 V 2vr 



(p-3p). (13) 



A plot of the quantity is shown in Fig. 2. Notice the linear growth. 

Asymptotically (A ^> A c ) this linear rise can be shown to hold analytically. To 
this end let us consider the function 

) = ^3p = ^3p = 1 + ?i ^3p M 
For a«l (i.e. A ^> A c ) we have [23] : 

- 71-2 1 2 - 71-2 1 2 

Pi o = 6 a + •••, p 19 = 6 -a + •••. (15) 

n,2 30 4 » fi,2 90 4 v ; 

As a consequence, the function h = h(X) becomes: 

ft W = 1 + jW + "' = 1 + T^i + "' = | + "-. ("I 

where Eq. (I12p has been used. Thus h(X) is approximately constant for a < 1. The 
numerical behavior for the function h = h(X) is plotted in the right panel of Fig.[IJ 
Notice that the asymptotic value h(X) = § is practically reached for A > 5 A c . 

By virtue of Eq. (1141) we have asymptotically (corresponding to the gray curve 
in the left panel of Fig[T]) : 

p - 3p = 6p gs = 24vrA 3 T . (17) 
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It is interesting that 9^ splits as = 6p gs = Ap gs + 2p gs where the first summand 
is the direct contribution of the ground state while the second summand arises from 
the massive modes a 1 ^ 2 . That is, the mass m = m(T) behaves in such a way that 
fluctuations generate a linear contribution to 9^ at high T. 

Let us now compare our prediction for the slope with lattice results. To this end 
we relate the Yang-Mills scale A to the critical temperature T c as 

A = ^£^0.45T C , (18) 

where A c = 13.89 has been used. Then we have asymptotically 

1 92 7T 4 

m = P-Sp= — — T 3 T ~ 7T C 3 T . (19) 

In the work of pQ a value of T c = 0.290 GeV was used. Thus Eq. (TT9~j) asymptotically 
predicts 9^ ~ 1.7 T GeV 3 . Notice that according to our approach the slope increases 
from ~ 1.4 GeV 3 to the asymptotic value 1.7 GeV 3 for 1.4 T c < T < 5T C . Thus 
our prediction for the slope is in agreement with that of Miller [1] whose lattice 
simulation yields ~ 1.5 GeV 3 (read off from his Fig. 1). In Ref. [3] a slope of ~ 
0.2 GeV 3 was obtained for the SU(2) case which also is in qualitative agreement 
with our result. Notice that the lattice simulations of [21] observes a linear behavior 
of 9^, too. 

In Fig. 2 we report the (scaled) energy density p, pressure p, and entropy density 
s as functions of temperature. A detailed comparison of these results with those 
obtained on the lattice (for both differential and integral method and for the SU(3) 
case) is carried out in [12J. While there is good agreement for the entropy density 
s (for which the ground state is not relevant) the pressure p becomes negative close 
to the phase boundary in virtue of the negative contribution of the ground state. 
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While this is in qualitative agreement with the result of the differential method, it is 
not with that of the integral method. This fact, as emphaiszed in the Introduction, 
represents a motivation to perform the present work where lattice results safely 
above the phase boundary, such as the behavior of W at high T, are compared with 
the effective approach under study. 



4 The SU(3) case 

In the SU(3) case the behavior of 9^ is qualitatively similar to the SU(2) case. We 
only report on some relevant formulas and briefly discuss their consequences. The 
modulus of the scalar field <fi is exactly the same. As shown in [12] out of the eight 
coarse-grained gauge modes four acquire a mass mi = e |0| (contributing to p and 
p by p\ and pi), two a mass m 2 = 2e \<p\ (p 2 an d P2), and two stay massless (p 3 and 
p 3 ). Explicitly, we have (a = rrii/T): 

A 6 f°° , xWx 2 + a 2 , N 

^ 3 = 4 30' * = 1? J dX e V ^-V (20) 

3 /"» x*y/x* + (2a)* _ 2(2tt) 4 
P2 = — / dx — , p gs = 3 , 21 



2 e poo 

p 3 = 4-, Pi = -- 2 J Q dxx 2 \n(l-e~^^y (22) 

p 2 = -- 2 I™ dxx 2 \n(l-e-^ 2+ ^ 2 ) , P 3S = -p 9S . (23) 
*/ 

Eq. ([H]) assumes the following form: 

1 = ~Wf + a ) {aD{a) + 2ajD(2a)) • (24) 

The asymptotic solution to Eq. (1241) reads a(A) = -^7r 2 A~ 3 / 2 , and the effective cou- 
pling reaches a plateau value of e = ^7r. The value for A c is A c = 9.475 [12] where 
the coupling and all masses diverge. The function h(X), defined in (1141) . has an 
asymptotic value of 3/2 like in the SU(2) case. Thus the asymptotically linear be- 
havior 9 up = p — 3p = 24nA 3 T holds. Relating A to the critical temperature T c as 
in Eq. ffT8l) . we take into account the SU(3)-value A c = 9.475. This yields: 

1927T 4 

0^ = P ~ 3p = ^^T C 3 T ~ 21.9 T 3 T . (25) 

As in Ref. [I] we set T c = 0.264 GeV and thus find 9^ ~ 0.4 T GeV 3 for T > T c . In 
[TJ a somewhat smaller slope of ~ 0.2 GeV 3 is extracted. In [3] a slope of 0.63 GeV 3 
was estimated. 
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5 Conclusions 



In this paper we studied the high-temperature behavior of the trace of the stress- 
energy tensor 9^ = p — 3p. This quantity measures the breaking of scale invariance 
and vanishes for a relativist gas of free gluons: recent lattice works [H[2] show a linear 
rise for high T (see however also [5J, where the rise is quadratic). This means that 
strong interactions still play a decisive role even in the high temperature domain. 

Theoretical studies should therefore be concerned with the predicted behavior of 
Oufi at high temperatures. We have shown that the linear growth of 9^ for high T 
is an analytic prediction of the approach to Yang-Mills thermodynamics [12j . which 
involves a thermal ground state emerging as an average over instantons at finite T. 
That is, even in the high-T limit the dynamic of these field configurations is impor- 
tant and can be detected by lattice simulations. In both cases, SU(2) and SU(3), we 
evaluated the slope of 9^ and the results 9^ ~ 7T C 3 T and 9^ ~ 21.9 T C 3 T , where T c 
is the deconfining temperature, are found for high T, respectively. We compared to 
the lattice results of Miller [TJ [2] and to the theoretical work of Zwanziger [3] finding 
good agreement. In the latter a momentum-dependent, universal modification of 
the dispersion relation for propagating, fundamental gluon fields, motivated by the 
reduction of the physical state space a la Gribov, is introduced. Our approach and 
[3] share the property of a nonpertrubative gluon mass leading to the linear growth 
of 9^. However, while in our case the mass is a temperature-dependent coarse 
grained quantity and the ground state itsself contributes to 9^, this is not the case 
in [3] . A more detailed comparison would be surely interesting but goes beyond the 
scope of the present work. 

We believe that future studies on the high T behavior of 9^ are therefore very 
important, both numerically and analytically, in order to confirm the linear growth 
and to relate it to nonperturbative aspects of Yang-Mills theories. 
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